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We investigate the possible role of linearly polarized gluons in Higgs production from 
unpolarized pp collisions. The transverse momentum distribution of the produced Higgs 
boson is found to exhibit a modulation with respect to the naive, unpolarized expectation, 
with the sign depending on the parity of the Higgs boson. The transverse momentum 
distribution of a scalar Higgs will, therefore, have a shape clearly different from a pseudo- 
scalar Higgs. We suggest that this effect can be used to determine the parity of the Higgs 
at the LHC, without the need to use challenging angular distributions. 

Introduction 

After a discovery of a new scalar particle at the LHC, the next task at hand is to determine its 
coupling to other particles. Not only the size, but also the type of coupling to fermions has to 
be determined, being either the P even <J"I' or the P odd ^75^. Relatively few suggestions to 
this end have been put forward for the LHC, e.g., using Higgs -I- 2 jet production jj or r pair 
decays [5] . We claim that the difference between a scalar and pseudoscalar coupling might also 
be visible in the transverse momentum distribution of the scalar particle. 

The Higgs transverse momentum distribution has been calculated in the framework of 
collinear factorization at Next-to-Next-to-Leading Logarithmic (NNLL) accuracy for small q T , 
matched to Next-to-Leading Order (NLO) accuracy for large q T [3111]. It was noted [H[S] that 
in NLO 77 continuum production there are "gluon spin-flip contributions" in the gg induced 
channel, which should be described by a "spin- flip distribution" V', that can, in principle, be 
as large as the unpolarized distribution. It was also noted [6] that in the NNLO radiative cor- 
rections to the Higgs boson q T cross section "gluon spin correlations" become important, which 
cause the standard Drell-Yan transverse momentum resummation to fail for the gluon-gluon 
fusion process. We claim that the fact that these gluon polarization effects are only observed at 
NNLO in Higgs production, is due to the use of the collinear factorization framework in which 
the polarized gluons have to be generated from the unpolarized distribution by gluon radiation. 
Within the framework of Transverse Momentum Dependent (TMD) factorization, the effect 
of polarized gluons is already present at tree-level and described by a non-perturbative input 
function h^ 9 . Although dependent on its size, the effects of polarized gluons are, in principle, 
large and modify the q T distribution of a scalar and pseudoscalar Higgs in a distinct way. 
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TMD factorization 

The contribution of gluon fusion to Higgs production in the TMD framework reads [7] in leading 
order in q T /rriH, 

C\(7 f / ~ \ / ~ \ * Pb — x hPb 

^ocj d 2 p T d 2 k T 5 2 (p T + k r - q T ) $£> a , p T )$ p /(x b , k r ) (M^j {M v °) p , (1) 

in which the momentum fractions are given by x a ^ — q 2 /{2P a ^ ■ q) and <I> g is the gluon 
correlator 0(5], 
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2x Y T J1 ^'^ T ' \ M 2 " T 2M 2 

with p 2 T = —p%, <7r" — 5 Miy ~ P^n v jP ■ n — n^P" / P ■ n, and M the proton mass. The function 
fi(x,Pr) represents the unpolarized gluon distribution and h-^ 9 (x,p%) represents the distribu- 
tion of linearly polarized gluons. 

Size of the linearly polarized gluon distribution 

As a first step, to study the effects of linearly polarized gluons, we follow a standard approach 
for TMDs in the literature and assume a simple Gaussian dependence of the gluon TMDs on 
transverse momentum: 

/f(s,p£) = ^exp(-if,Y (3) 



^(Pt) V (p 

where G(x) is the collinear gluon distribution, G(x) = J d 2 p T ff (x, p T ). The width, (p 2 ) 7 
depends on the energy scale, Q, and should be experimentally determined. We will estimate 
(p 2 ) = 7GeV 2 , at Q = mn — 120 GeV, in rough agreement with the Gaussian fit to ]^(x,p 2 T ) 
evolved to Q = M z of Rcf. [TO]. 

No experimental data on hf 9 is available, but a positivity bound has been derived in Ref. [8J: 
\hi 9 (x,p%.)\ < f 9 (x,p 2 ). We will use a Gaussian Ansatz for h^ 9 , with a width of r(p\), 

. ±g , 2 , M 2 G(x) 2e(l - r) / p\ \ 

with a normalization such that it satisfies the bound for all p T . 
Higgs transverse momentum distribution 

Using the TMD factorization expression in Eq. (UJ), the parameterization of the gluon corre- 
lator in Eq. @ and the Ansatz for the TMD distribution in Eqs. and (UJ), the transverse 
momentum distribution for a scalar and pseudoscalar Higgs can be written as 
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where ± stands for a scalar/pseudoscalar and R(q T ) = C[wHhf 9 hf 9 ]/C[ffff], with 



d 2 p T / d 2 k T 5 2 {p T + k T - q T )w(p T , k T ) f(x a , p 2 ) f(x b , k 2 ) 



C[wff] 



and wh — [(p T ■ k T ) 2 — ip 2 k 2 ]/2A/ 4 . With our Ansatz for the distribution functions, 
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exp 



1 



2(P 2 ) 



(6) 



(7) 



r{p 2 T ) 8r 2 (p2} 2 , 

The transverse momentum distribution in Eq. ([5]) for a scalar and pseudoscalar Higgs is plotted 
in Figure Q] for r — 2/3 and r = 1/3. As long as h^ g is not measured, the absolute size of the 
effect is unknown, but it will always be such that a scalar has enhancement at low q T , suppression 
at moderate q T , followed again by enhancement at high q T , whereas for the pseudoscalar this 
is reversed. Higher order perturbative corrections will modify the exact form and width of our 
tree-level q T distribution, as well as the size of the modulation, but this qualitative behavior is 
expected not to change. 
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Figure 1: Transverse momentum distribution of the Higgs, using the parameterization of h 1 
in Eq. Q with r = 2/3 (left) and r = 1/3 (right). The naive curve is the prediction for both 
scalar and pseudoscalar in case h 1 9 = 0. 



Two photon decay channel 

When the Higgs decays to, e.g., two photons, there will be irreducible background due to 
gg — > quark box — > 77, which was recently investigated in the framework of TMD factorization 
|11) . Including this background, we come to the conclusion [7], that the transverse momentum 
distribution of the photon pair has the same form as Eq. ([5]), but with a Q and Collins-Soper 
angle, 9, dependent size, i.e., 



1 da 
a dq 2 



1 



F 2 (Q,0) 



R(q T ) 



(8) 



Fi(Q,e) 

The ratio F2/F1 is plotted in the left graph of Fig. [5] At the Higgs mass we reproduce Eq. 
([5]), i.e. F2/FI — > ±1, for a scalar/pseudoscalar, but away from the pole, the background 
quickly dominates. To mimic a finite detector resolution in the determination of Q, we also 
plot the ratio F2 / F% in which both numerator and denominator are separately weighted with a 
Gaussian distribution. From the graph we see that the continuum background reduces the effect 
to approximately 30% or 20% of the maximal size with a 0.5 or 1 GeV resolution, respectively. 
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Figure 2: The ratio F2/F1 in Eq. ([8]) plotted as function of Q for 9 = ir/2 assuming a 120 GeV 
Higgs (left) and the same curves including a detector resolution of 0.5 and 1 GeV (right). 

Conclusions 

The effect of gluon polarization in pp collisions is such that scalar and pseudoscalar particles, 
produced through gluon fusion, will have different transverse momentum distributions. Al- 
though the absolute size of the effect cannot be estimated without experimental input on h^ 9 , 
the qualitative features are such that this effect can, in principle, be used to distinguish scalar 
from pseudoscalar particles. In the two photon decay channel of the Higgs, the continuum 
background partially washes out the difference between scalar and pseudoscalar. Other decay 
channels are currently being investigated. 

This work is part of the research program of the "Stichting voor Fundamenteel Onder- 
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